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Abstract 

We analyze a method to produce pairs of non independent Poisson pro¬ 
cesses M{t),N{t) from positively correlated, self-decomposable, exponential 
renewals. In particular the present paper provides the family of copulas pair¬ 
ing the renewals, along with the closed form for the joint distribution Pm,n{s, t) 
of the pair [M{s),N{t)), an outcome which turns out to be instrumental to 
produce explicit algorithms for applications in finance and queuing theory. 
We finally discuss the cross-correlation properties of the two processes and 
the relative timing of their jumps 


1 Introduction 

Recent studies have shown that the spot dynamics of commodity markets displays 
mean reversion, seasonality and jumps [1], and some methodologies have been pro¬ 
posed to take dependency into account based on correlation and co-integration [2]. 
However, these approaches can become mathematically cumbersome and non-treatable 
when leaving the Gaussian-Itd world. In this context it has been indeed recently 
proposed [3] to consider 2-dimensional jump diffusion processes with a 2-dimensional 
Gaussian and a 2-dimensional compound Poisson component, and, as also suggested 
in different circumstances [4], we show here that a revealing approach to model the 
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dependency of the 2-dimensional Poisson processes can be supplied on the ground of 
the self-decomposability of the exponential random variables used for its construc¬ 
tion. 

The present paper is in particular devoted to hud both the copula function pair¬ 
ing our correlated renewals, and an explicit form for the joint distribution Pm.,n{s, t) of 
our pair of correlated Poisson processes M{s), N{t): this will prove to be instrumen¬ 
tal to produce the efficient algorithms that can be used in hnancial applications [3]. 
If indeed the pairs of correlated, exponential random variables (rt^’s) (X^, Y^) - used 
to produce the renewals in our processes - are interpreted as random waiting times 
with random delays, the proposed model can help describing their co-movement and 
can answer some common questions arising in the hnancial context: 

• Once a hnancial institution defaults how long should one wait for a dependent 
institution to default too? 

• A market receives a news interpreted as a shock: how long should one wait to 
see the propagation of that shock onto a dependent market? 

• What is the impact of the correlations among the shocks for diherent insurance 
companies on a fair assessment of the risk of losses [5] ? 

It is worth noticing, moreover, that we achieve our aim of producing a 2-dimensional 
Poisson process with dependent marginals without resorting to an a priori copula 
(distributional) approach: the dependence among arrival times will indeed be made 
explicit in terms of combinations of rv’s, and we only recover and discuss the corre¬ 
sponding copula functions as an outcome of this model. As a consequence, because 
of this P-a.s. relationship between the random times, the two Poisson processes can 
be seen as linked with a form of co-integration between their jumps. Similar models - 
albeit rather less sophisticated - were also used in order to model a multi-component 
reliability system [4], while the so-called Common Poisson Shock Models [6] are in 
fact quite different from that presented here 

The main practical consequence of our results is then that the price and the 
Greeks of the spread options considered in the applications [3] can be calculated in 
closed form using either the Margrabe formula (if the strike is zero), or some well 
known approximation [7]. In any case our model entails explicit algorithms for the 
simulation of correlated Poisson processes, and can be used in the Monte Carlo simu¬ 
lations. An extension to the multi-dimensional case, as well as to different dynamics 
other than Poisson, will be considered in future studies; but, under the assumption 
that only two underlyings have jump component, the price and the Greeks of spread 
options can be obtained even now by the moment-matching methodology recently 
proposed in [8] 

The paper is organized as follows: in the Section 2 we hrst show how (hitherto 
positively) correlated exponential rv^s can be deduced from the self-decomposability 
of their laws; then in the Section 3 we briefly discuss the copula functions produced 
by this model. By using pairs of these exponential rv^s as correlated renewals, in the 
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Section 4 we subsequently produce a 2-diniensional Poisson process with correlated 
components, and in the Section 5 we explicitly deduce their joint distributions. 
Finally in the Section 6 the cross-correlation properties of the Poisson processes are 
briefly analyzed, and the relative timing of their jumps is used to shed new light 
on the dependence mechanism of a model allowing for the possibility of a delayed 
propagation of correlated shocks. We conclude by pointing out hrst that we would 
also be able to produce other correlated rv^s (Erlang, Gamma, EPT...) by making 
use, once more, of their self-decomposability; and then that the results of this paper 
should also be extended to negatively correlated renewals, a possibility - not open 
to other procedures - that will be postponed to future inquiries. Lengthy proofs are 
conhned in the Appendices, together with a few details about the notation adopted 
throughout the paper 

2 Correlation from self-decomposability 

2.1 Joint distributions 

A law with density {pdf) f{x) and characteristic function {ehf) (p{u) is said to be 
self-deeomposable {sd) [9, 10] when for every 0 < a < 1 we can End another law 
with pdf ga{x) and ehf Xa{u) such that 

cp{u) = ^p{au)xa{u) 

This is a well known family of laws with many relevant properties. We will also say 
that a random variable {rv) X with pdf f{x) and ehf ip{u) is sd when its law is sd: 
looking at the dehnition this means that for every 0 < a < 1 we can always find 
two independent rv^s Y (with the same law of X) and Za (with pdf ga{x) and ehf 
Xa{u)) such that 

X = aY + Za 

We can look at this, however, also from a different perspective: if Y is sd, and to 
the extent that, for 0 < a < 1, an independent Za with the suitable law is known, 
we can define a third rv 


X = oY + Za P-a.s. 


being sure that it will have the same law as Y. In the following we will mainly adopt 
this second standpoint 

We turn now, for later convenience, to give the joint laws of the triplet {X, Y, Za): 
for the ehf 'ip{u,v,w) we easily find from the independence of Y and Za that 


'lp{u, V, w) 


^^i{uX+vY-\-wZ 


a) 


ip{au + v)xa{u + w) = ip{au + v) 


ip{u + w) 
(p{a{u + w)) 
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while the marginal, joint ch/’s of the pairs {X, Y) and {X, Za) respectively are 


(t){u,v) 
u{u, w) 


V, 0) = (p{au + V 
0, w) = ip{au 


Pju) 
ip{au) 
ip{u + w) 


(p{a{u + w)) 


As for the pdf K{x,y,z) on the other hand , by taking s = au + v, we have (we 
neglect the integration limits whenever they extend to all R) 




(27r)3 

1 

1 


duJdvjd.e-^--^y--^^iau + v)UuP.) 

^du jdw + ^)^ / + v) dv 


= fiy)7YXJ du + 


(27r)2 . 

and again with s = u + w, t = w 

1 


K(x,y,z) = f(;„)—jds y<i«e-‘K-‘‘»><-'>+-'‘lxa(s) 


27r 


27r 


= /(d) ga{,x - ay) S[z -{x- ay)] 


so that the marginal, joint pd/’s of (X, Y) and (X, Za) will respectively be 

h{x,y) = f{y)ga{x-ay) i(x,z) = - f(- — ga{z) (1) 

a \ a J 

Finally the joint cumulative distribution function (cd/) of (X, X) is 

H{x,y)= f f{y')Ga{x - ay')dy' Ga{z) = f ga{z') dz' 

J —oo J —oo 

where Ga{z) is the cdf of Za- The particular form of H{x,y) will be instrumental 
in hnding the copula functions [11] eventually pairing X and Y 

We can hnally also calculate the correlation coefficients rxr and rxzZ if pat 
E [X] = E [Y] = fi and W [X] = V [Y] = from the Y, Za independence we have 

E [XX] = E [(aX + Za)Y] = aa^ + 


and hence txy = a. In a similar vein, to calculate rxZa we hrst remark that 
V [X] = a'^V [X] + V [Za], namely V [Za] = (1 — af)a‘^, and then from 

E [XZa] = E [{aY + Za)Za] = (1 - a^)a^ + (1 - a)y^ 


we hnally hnd rxZa = 1 — 
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2.2 An example: the exponential laws ^i(A) 

It is well known that the exponential laws Si (A) with pdf and chf (see Appendix A 
for the notations adopted from now on) 

A/i(Aa;) = \e~^^d{x) Pi{u/\) = ^ . 

A — lU 


are a typical example of sd laws [9], and in this case we can explicitly give the law 
of Za- we have indeed 


Xa{u) = 


(pi{u/\) A — iau 
ipi{au/\) \ — iu 


— O "p (1 — o) 


A 


A 


lU 


= a+{1 - a)(pi{u/\) (2) 


which (for 0 < a < 1) is a mixture of a law (5o degenerate in 0, and an exponential 
Si (A), namely 

Za ~ ado + (1 — a)Si(A) 
so that its pdf and cdf respectively are 


Qaiz) = a6{z) + (1 — a)Ae 
Ga{z) = [a + {l-a){l-e~^^)]d{z) 


It is also easy to prove on the other hand that this coincides with the law of the 
product of two other independent rv^s: an exponential Z ~ Si (A), and a Bernoulli 
B(l) 03(1,1 — a) with a = P {-B(l) = 0}, so that we can always write 


Za = B{l)Z 

In short, given two exponential rv^s Y ~ Si(A) and Z ~ Si(A), and a Bernoulli 
5(1) r\-/ iB(l, 1 — a), all three mutually independent, the rv X dehned as 


X = aY + B{l)Z (3) 

is again an exponential Si (A). From (1) we also hnd that the joint pdf of X, Y is 
h{x,y) = Xe~^'^d{y) \a5{x — ay) + (1 — — ay)] 

and hence its joint cdf is 

H{x, y) = r Xe-^y'd{y') [a + (1 - a)(l - d{x - ay') dy' 

= 1) (s A ^) [ (l - (l - )) 


Of course this is far from the only possible joint law with exponential marginals (see 
also Section 3.1), but it is noticeable because it traces its origins back to a model 
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Figure 1: The pairs (Xk, Yk) with correlation 0.01 and 0.99 


of self-decomposability of the exponentials. As for the correlations among X, Y and 
Z, we already know that txy = For we hrst hnd that 

E [XZ] = E [{aY + Za)Z] = ^ 


and then that 


rxz = 1 - a = 1 - rxY 


so that for our three exponentials in (3) we eventually have 


rxY + rxz = 1 rxY = a ryz = 0 


2.3 Positively correlated exponential rv’s 

It is apparent now from the discussion in the previous section that the self-decom- 
posability of the exponential laws £i(A) can be turned into a simple procedure to 
generate identically distributed and correlated rv’s: given Y ~ '£i(A), in order to 
produce another X ~ Ci(A) with correlation 0 < a < 1, h would be enough to take 
Z ~ Si(A) and B{1) ~ ®(1,1 — a) independent from Y and dehne X as in (3). In 
other words X will be nothing else than the exponential Y down a-rescaled, plus 
another independent exponential Z randomly intermittent with frequency 1 — a. The 
self-decomposability of the exponential laws ensures then that, for every 0 < a < 1, 
also X marginally is an £i(A). Remark that we would not have the same result 
by taking more naive combinations of Y and Z. Consider for instance the sum 
oY + {1 — a)Z of our two independent, exponential rt>’s: in this case, since oY ~ 
(A) and (1 —a)Z ~ fhe law of aY + {l — a)Z would be (A) (i^); 

which is neither an exponential Si (A), nor even in general an Erlang S 2 because of 
the difference between the two parameters 

The proposed procedure can now be adapted to generate a sequence of indepen¬ 
dent pairs of exponential rv’s, with correlated components, (Xfc,Yfc), k = 1,2,... 
(or, if we prefer, Xk,Zk) that will act in the subsequent sections as renewals for 
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a two-dimensional point (Poisson) process: take indeed 0 < a < 1, produce two 
independent id exponentials Y, Z and another independent Bernoulli -B(l), then de- 
hne X = aY + B{1)Z and take the pair X,Y. By independently replicating this 
procedure we will get a sequence of iid two-dimensional pairs {Xk,Yk) that will be 
used later to generate a two-dimensional Poisson process with correlated renewals. 
Of course - not surprisingly - the case of uncorrelated pairs of renewals {Xk,Yk), 
and hence of independent Poisson processes, is retrieved from our model in the limit 
a —)■ 0, because in this case we just have X = Z which is by dehnition independent 
from Y. On the other hand it is also apparent that in the opposite limit a —)■ 1 
(namely when txy goes to 1) we tend to have X = Y, P-a.s. so that the time pairs 
will fall precisely on the diagonal of the two-dimensional time, and the two Poisson 
processes will simply P-a.s. coincide. To see it from another standpoint we could 
look to some simulation of the pairs (X^, Yk): for small correlations a the scatter-plot 
of our pairs {Xk, Yk) tends to evenly spread out within the hrst quadrant without 
any apparent hint to some forme of dependence; on the other hand for a near to 1 
the points tend to cluster together along the diagonal, as can be seen in the Figure 1 

3 Copulas for bivariate exponentials 

3.1 A family of copula functions 

From the discussion in the Section 2.2 we know that the pair X, Y of correlated rv^s 
deduced from their self-decomposability has the joint cdf 

H{x,y) = d(yA'^'^ (l - (l - )) (4) 

with the exponential marginal cd/’s (the notation is here slightly simplihed) 

F{x) = d{x) (1 - e-^^) G{y) = d{y) (l - (5) 

To hnd out the copula function C{u, v) [11] pairing X, Y we then first remark that 

g-Ax _ ^ ^-Xy ^ ^ ^-aXy _ _ Q^y^^a 

while 

0<-<y ^ ^ [1 - G{y)f < 1 - F{x) 

ct 

0<y<- ^ ^ [1 - G{y)f' > 1 - F{x) 

(X 

and then that our joint cdf (4) takes the form 


H{x,y) = F{x) 

H(x,y) = F(x) - [1 - G(s)| ([l - 


for [1 — G(|/)]“ < 1 — F{x) 
for [1 — G(|/)]“ > 1 — F(a;) 
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which can also be conveniently summarized as 


H(x,y) = F(x)-[l-G(y)] 


l-F(x) \ 


As a consequence we get the following family of copula functions 


Ca(u,v) = U - (1 -V) 



1 — U 

(1 — n)“ 


[(1 — n )“ — (1 — u )]^ 
(1 — 


which for 0 < a < 1 runs between two extremal copulas 


( 6 ) 


Co{u,v) = uv independent marginals 

Ci{u,v) = uAv fully positively correlated marginals 

It is easy to see that Ci{u,v) also coincides with the Frechet-Hoffding upper bound 
C{u,v) for copulas (see Section 3.3) 


3.2 Bivariate exponential distributions 

Several examples - all different from (6) - of bivariate distributions with exponential 
marginals Si(A) and can be found in the literature [11, 12], First we hnd the 

Gumbel bivariate exponential distribution [11, 13] with 0 < a < 1 and 

H{x, y) = d{x)d{y) (l - 

h{x, y) = d{x)d{y) [a(Aa; + yy + aXy xy) + 1 — a] 

Gaiu.v) = M + + 

It is apparent that Gq{u,v) = uv gives the independent exponentials, while 
Gi{u,v) = U + V — 1 + {1 — U){1 — V)e~^ In(l--u) ln(l-t>) 


does not seem to correspond to some notable copula. Then there is the Marshall- 
Olkin bivariate exponential distribution [11, 14] with 0 < a, 6 < 1 and 


H{x,y) 

h{x,y) 


f d(x)d(y)(l - e-^^y-^(l - if (1 - > (1 - 

{d(x)d(y)(l - e-^^)(l - if (1 - < (1 - 

J (i_eAL)a 'd(a;)Ae“'^* d{y)ye~^^'^ if (1 — > (1 — 

1 {i-y-\y)b 'd{x)Xe~^^ 'd{y)ye~^y if (1 — < (1 — 


Ga,b{u,v) = {u^ “n) A {uv^ y 




when 

when 


In this case Gofi{u,v) = uv again is the independent copula, while Gi^i{u,v) = uAv 
is the Frechet-Hoffding upper bound G{u,v) (see Section 3.3): apart from these 



N Cufaro Petroni and P Sabino: Correlated Poisson processes 


9 


extremal values, however, also this Marshall-Olkin copula differs from (6). A third 
family of copulas can finally be traced back to the Raftery bivariate exponential 
distribution [11, 15]: in this case the copula functions are 


Ca{u, v) =uAv + —^{uv )' 
2 — a 


(uV V 


a-4 


and correspond to the case of correlated exponential rv^s X,Y which are produced 
by three independent exponential rv^s U, V and Z according to the definitions 


X = aU + B{l)Z 


Y = aV + B{l)Z 


Here, at variance with our model based on self-decomposability, the correlation is 
apparently produced by the presence of the same exponential rv Z in both the right- 
hand sides of the definitions. In short, it results from these examples that our family 
of copulas (6) seems not to have been used in advance to couple pairs of marginal 
exponentials 


3.3 Frechet-Hoffding bounds 

It is well known known [11] that every copula function C{u,v) falls between the 
Frechet-Hoffding bounds 

C_{u, v) = {u + v — !)■'■ < C{u,v) < u Av = C{u,v) 


and we have also found in the Section 3.1 that the copula Ci{u,v) for our fully 
correlated {rxY = 1) exponential marginals coincides with the Frechet-Hoffding 
upper bound. By keeping in mind a possible generalization of our model to the 
case of negatively correlated exponentials, we will briefly recall in this section a few 
general features of the joint cd/’s H{x, y) = C{F{x), Giy)) produced by the pairing 
of two given cd/’s F{x) and G{x) by means of the Frechet-Hoffding lower and upper 
bounds 

Let us suppose for simplicity that F{x) and G{x) are continuous and strictly 
increasing functions so that the inverse functions exist, and consider first the lower 
bound copula C_{u, v) = {u + v — I)'*'; in that case the condition F{x) -|- G{y) > 1 is 
equivalent to both the inequalities 

X > (d{y) = F-\l - G{y)) y > a{x) = G-\l - F{x)) 


and hence from H{x,y) = {F{x) + G{y) — 1)+ we first have 


d^H{x,y) 

dyH{x,y) 


fix) 

if 

y > aix) 

0 

if 

y < oi{x) 

aiy) 

if 

X > (3{y) 

0 

if 

X < I3{y) 
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where f(x) and g(y) are the corresponding marginal pd/’s, and then the joint pdf is 

h{x, y) = d^dyH{x, y) = f{x)5{y - a{x)) = g{y)5{x - f3{y)) (7) 

As a consequence we can say that the joint laws produced by the copula C{u, v) 
describe pairs of coupled rv^s X,Y satisfying P-a.s. the functional relations 

X = I3{Y) = F-^{l-G{Y)) Y = a{X)=G-\l-F{X)) (8) 

A formally identical result can be proved for the case of the upper bound 2-copula 
G{u,v) = u A V but for the fact that now the functions a{x) and /3{y) must be 
redehned as 

a{x) = G-\F{x)) m = F-\G{y)) 

In the case of the lower bound copula C{u, v) it is interesting to remark now 
that when for instance F{x) and G{y) are Gaussian cd/’s the functions a{x) and 
I3{y) are linear with negative proportionality coefficients, so that the pair X,Y is 
perfectly anti-correlated with txy = —1- The same happens also in the case of a 
pair of Student laws of the same order. This is true indeed for every other pair of 
marginal laws of the same type and with support coincident with R. On the other 
hand when the marginals either are not of the same type, or have an unbounded 
support strictly contained in R (as happens for exponential laws), they apparently 
can not reciprocally be in a linear relation with negative proportionality coefficient, 
and hence can not be totally linearly anti-correlated. In this case it can still be 
proved by means of Hoffding’s Lemma (see [11] p. 190) that the minimal correlation 
is reached by means of the lower bound copula C, but now a{x) and (diy) can 
no longer be linear functions, and rxv will be strictly larger than —1. By taking 
indeed the Frechet-Hoffding lower bound Cfu, v) = {u + v — I)’*' as the copula for our 
exponentials (5) we would hud the pdf (7) and the functional relations (8) where 
now 

aix) = -j In (1 - e-"") ^ = — In (l - 

and a short calculation would then show that in this case 

TT^ 

rxY = 1-~ —0.645 

6 

so that this minimal anti-correlation allowed for exponential rv’s would in any case 
be larger than —1. It could in fact be proved in general (see [11] p. 30-32) that, 
when X and Y are continuous, Y is almost surely a decreasing function of X if and 
only if the copula of X and Y is C_. Random variables with copula C_ are often 
called countermonotonic. We postpone to a subsequent enquiry a detailed study of 
negatively correlated exponentials 
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Figure 2: Sample pairs of the two-dimensional point process (T„, Sn) with correlation 
rxY = 0.01: on the left the points are compared with the average trend ^); on 
the right they are instead plotted after centering around these averages 


4 Correlated Poisson processes 

Following the discussion of Section 2 it is easy now to produce a sequence of rt>’s by 
independently iterating the dehnition (3) 

Xk = aYk P k = l,2,... (9) 

in such a way that for every k\ are (Si(A), i?fc(l) are 53(1,1 — a), and 

Yk, Zk, Bk{l) are mutually independent. The pairs {Xk,Yk) instead will be a-corre- 
lated for every k. Add moreover Xq = Yq = Zq = 0, P-a.s. to the list, and take 
then the point processes for n = 0,1, 2,... 

n \ ^ ^ 

= = Rn = J2Zk ( 10 ) 

k=0 ^ fc=0 k=0 

Since the Xk ~ '£i(A) are iid rv’s we know that Tn ~ Sn(A) are distributed as 
Erlang (gamma) laws with pd/’s Xfn{Xx) and ch/’s ipn{u/X) (see Appendix A for 
notations) where it is understood that Tq ~ Co = ^o- In a similar way we can argue 
that Sn ~ Cn(/i) and ~ Cn(A). We will hnally denote with N{t) ~ ^(At) and 
M(t) ~ the correlated Poisson processes associated respectively to Tn and Sn 

In order to get a hrst look to these processes we generate n = 1 000 pairs (X^, W) 
with the associated two dimensional point process (T„, Sn), and then we simulate the 
corresponding Poisson processes N(t) and M(t). The pairs (T„, Sn) are hrst plotted 
along with their average time increases (f, ^), and then after centering around these 
averages, namely as 

n ^ n 

Tn-- Sn-- n = 1 , 2 ,..., 1000 

A pi 
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Figure 3: Sample pairs of the two-dimensional point process (T„, Sn) with correlation 
rxY = 0.99: on the left the points are compared with the average trend ^); on 
the right they are instead plotted after centering around these averages 


In this second rendering the random behavior is magnihed by consistently reducing 
the plot scale to a suitable size. In the same way for the Poisson processes we hrst 
show samples of the pair and then that of their compensated versions 

N{t) = N{t) — Xt and M{t) = M{t) — fit 

In the Figure 2 we plotted the two-dimensional point process (T„, Sn) with 
X = fi = 1 and a = txy = 0.01: since the correlation among the renewals is 
negligible the right-hand plots (centered around the averages) apparently show a 
random behavior. In the Figure 3 instead we took a = txy = 0.99, namely we 
generated strongly and positively correlated renewals. In this second case, as it was 
to be expected, the centered time pairs fall into line among themselves. As for the 
Poisson processes themselves, in the Figure 4 the trajectories on the left hand side 
have a = 0.01 correlation and look fairly independent, after suitable compensation, 
on the right hand side. In the Figure 5 instead we took a correlation a = 0.99 
and the compensated trajectories are now almost superimposed. Remark as on the 
left-hand sides of these hgures both the Poisson processes and the time pairs appear 
to be quite near to one another, and to their averages because of a scale effect which 
is eliminated by compensation and centering in the corresponding right-hand sides 

Proposition 4.1. The rv’s 

n 

Cn = ^ Bk{l)Zk 
k=0 

turn out to he the sum of a (random) binomial number B{n) ~ IB(n, 1 — a) of iid 
exponentials Ci(A), and hence they follow an Erlang law with a random index B{n) 
(here B{0) = Oj, namely 

B(n) 

Cn = ^ Zk ~ €s(„)(A) 

k=0 
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Figure 4: On the left sample paths are shown of the two Poisson processes N(t) and 
M(t) with correlation rxr = 0.01; on the right we instead have the corresponding 
compensated Poisson processes N(t) and M(t) 


Proof: This is better seen from the point of view of the mixtures by remarking 
that, if cpi(u/A) is the chf of £i(A), we have from (2) (see also Appendix A) 


n 


n 






) = E = E 


.k=0 






iuBkil)Zkl 


k=0 


k=0 ^ ' 


k=0 


which, if (pk{u) = [(pi{u)Y are the chf of <£fc(l), eventually is a mixture of Erlang 
laws iSfc(A) with the binomial weights It is understood here that = 1, 

so that £o(A) = (5o and /o(a:) = 5{x) (see Appendix A) □ 


A straightforward consequence of the previous proposition (which apparently just 
amounts to acknowledge a subordination) is that now from 


n n 




A:=0 


k=0 k=0 


we will also have 


B{n) 

rp _ C A _ o I \ ^ rpr _ C D 

J-n ^ “T ^ ~r ^ ^ ^ \ ^B{n) 

k=0 


( 11 ) 


where RB{n) ~ ^B(n)(A) is the point process Rn with a random index B{n). It is 
worthwhile to notice that the previous results also substantiate the well known fact 
that the Erlang rv’s are self-decomposable too: the explicit knowledge of the C,n law 
allows indeed to construct pairs of dependent Erlang rv’s with correlation a 












N Cufaro Petroni and P Sabino: Correlated Poisson processes 


14 



Figure 5: On the left sample paths are shown of the two Poisson processes N(t) and 
M(t) with correlation rxr = 0.99; on the right we instead have the corresponding 
compensated Poisson processes N(t) and M(t) 


5 The joint distribution 

Our main task is now to explicitly calculate the joint distribution of our Poisson 
processes at arbitrary times s,t > 0 and n, m = 0 , 1, 2 ,... 

Pm,n(s, i) = P{M (s) = m, N{t) = n} 

= P {5'm < S < S'm+1, Tn<t < Tn+l} 

and to this effect we hrst remark (in a slightly simplihed notation) that 

Pm,n = P {M{s) > m, N{t) > n} — P {M{s) > m + 1, N{t) > n} 

—P {M{s) > m, N{t) > n + 1} + P {M{s) > m + 1, N{t) > n + 1} 

~ Qm,n Qm,n+1 4“ Q'm+l,n+l (1^) 


where 

gm,n(s,t) = P{M{s) > m, N{t) >n} = P {Sm < S,Tn< t} 
so that by taking 

Ar Xt Xt — aps 

w = — y = — ^ =-< y 

a a a 


from (11) we are reduced to calculate (see also Appendix A) 
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Qm,n P "I S'm — ^ Sn P RB{n) — ^ 

^ poo 

= '^l^ein) / dr\fe{\r) 


e=o 


I *S'm < -S, ^Sn + RB{n) <t Ri = r, B{n) = 


^ X] / drfi{\r)P < S'm < s, S'n < A 


e=o 

ry 


t — r 


y- z 


ay 

y-w 


= a dw hn{aw)P < Sm < -, Sn < 

Jo [ P P 

where Xhn{\x) is the (Erlang binomial mixture) pdf of RB{n) 

Proposition 5.1. For n,m = 0,1,2,... and p,T > 0 we have 
P {Sm< p, Sn<T} = Pm\/n{y{p P t)) 

(mVn) —1 

+ [Qn-ra'&{T - p) + Qra-n'&{p - r)] ^ n(mVn)-fc (h|P “ t| ) TTfc (/i(p A t)) 

k=mAn 


(13) 


with the notations adopted in the Appendix A for the Poisson laws 

Proof: See Appendix B for a detailed proof □ 

Of course in (13) we take in particular 

y — z y — w .t — r 

p = - = s T = - = A- 

p p ay 

It is apparent that this result will be instrumental to calculate hrst qm,n{s, t) in (13), 
and then the distributions Pm,n{s,t) 

Proposition 5.2. If ays > \t, then Pm,n{s,t) = 0 whenever m < n 

Proof: Since our renewals Xk,Yk,Zk are all non-negative rvA, the point processes 
are always non-decreasing 


Sm < Tn < Tn+i m,n = 0,1,2,... 


while 


r„ = f s„ + 5: z. > ((f s. 


B(n) 


ay ^ 


k=0 


Now, if M{s) = m and N{f) = n, we must have both Sm < -s < Sm+i and <t < 

T„+i. Suppose now 0 < m < n, namely m -|- 1 < n and Sm+i < Sn- then 


ay ay ^ q ^ ^ . 

^ >5 \ ^ iJm+l A ^ S k 
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which apparently contradicts the hypothesized inequality □ 

As a consequence when a/is > Xt we can always restrict our calculations to the 
case m > n > 0. We can now dually state our complete results about the joint 
distributions Pm,n{s, t) 


Proposition 5.3. Take for short 

Xt Xt — aas 

y = — >0 =-< y 

a a 

Then, when ays > Xt, namely z <b, we have 

{ 0 n > m > 0 

Qn,n{y,z) m = n>b (14) 

Qm,niy 1 z) Qm,n+l(yj z') 771 ^ n P 0 

where we defined 

fy m ^ 

Qm,n{y^^)=(^ dwhn{aw)'Yl'Km-k{w - z)r:k{y - w) m>n>0 (15) 


When instead ays < Xt, and hence 0 < z < y, we have 
Pm,n{,y-i z) 

{ Am,n(y, Z) - Arn^n+l(y, z) + z) - 5„,n-l(l/, z) n > 771 > 0 

An,nidJi z) An^n+l(yj z') p z') A z') 771 77 A 0 

Am,n{y^ - Am,n+l{y, z) + Cm,n{y, z) - Cm,n+l{y, z) 777 > 77 > 0 

where we defined 


Am,n{,y, z) = a dw hn{aw)nm{y - z) 


77 , 777 > 0 


Bm,n{y, z) = a dw hn+fiaw) '^TTk{z - w)TTmiy - z) 77 > 777 > 0 (18) 

k=0 

py ^ 

Cm,n{,y,z) = a dw hn{aw)'^Tl„^-k{w - z)Tlk{,y - w) 777 > 77 > 1 (19) 


while for m > n = 0 we always have Cmfi{y,z) = 0. Moreover both the results for 
z < 0, and for 0 < z < y connect with continuity in z = 0 in the sense that 

Pm,n(l/,0") = Pm,n(7/,0+) 777 , 77 > 0 

Proof: Take drst the case ays > Xt, namely z < 0, and recall that for the integra¬ 
tion variable in (13) it is 0 < tc < 7/. As a consequence, when the Proposition 5.1 
in used in (13), we will always have 

y-^ ^y- ^ (oo\ 

0 < r =-<- = s = p (20) 
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On the other hand, since the conditions of the Proposition 5.2 are met, we can also 
restrict ourselves to evaluate Pm,n(s,t) for 0 < n < m. Then, by considering sepa¬ 
rately the cases m = n > 0 and m > n > 0, from (13) and from the Proposition 5.1 
we hrst calculate gm,n, Qm+i,n, Qm,n+i and gm+i,n+i, and hnally (lengthy algebraic 
details can be found in the Appendix C) from (12) we hnd (32). 

When on the other hand aps < Xt (namely y > z > 0 and 0 < w < y) and we 
use Proposition 5.1 in (13), instead of (29) we hnd 


0 < r 


y-w 

y 


0 < p = s = 


y- z 
y 


( 21 ) 


SO that p and r can now be in an order whatsoever. As a consequence Proposition 5.2 
does not hold, and we must consider all the possible orderings of m, n. Following 
then the same line of reasoning as before, and always taking separately the different 
n, m orderings, a tedious calculation (see Appendix C) gives hrst the g’s from (13), 
and eventually the p’s of our proposition from (12) 

We hnally show that the values of Prn,n{yiz) separately listed in the Proposi¬ 
tion 5.3 for z < 0 and z > 0 connect with continuity in = 0, in the sense that for 
every y > 0 

Pm,nijji 0 ) Pm,n{yi 0 ) 

For z < 0 (namely ays > Xt) the results are given in (32) and (15), so that for 
z 0~ and every m, n > 0, we simply have 


Pm,n(l/,0 ) 

Pn^nigii 0 ) 

Pm,n(|/,0") 


= 0 


= a dw hn{aw) 7io{w) 7in{y — w) 


n > m >0 

(22) 

n = m > 0 

(23) 

— w) m > n >0 

(24) 

k{w)7ik{y - w)l 



k=n 


On the other hand, when z > 0 (namely ays < Xt) we have (16), (17), (18) and (19), 
so that now appears also as an integration limit, and some care should be exercised 
for O’*". When indeed the integrand contains the distribution S{x), as in fact hap¬ 
pens in every hrst term of h„(a;) which is /3o{n)S{x) = a^5{x) (see also Appendix A), 
we have for every regular function ^{x) 


lim [ ^(x)6(x) dx = ^(0) 

2:4.0+ Jq 

As a consequence we have 

lim / dx^{x)hn{x) = 
Z4.0+ Jq 


ry 

lim / dx^(x)hn(x) = 

zl0+ 


ry 

lim / ^{x)5(x) dx = 0 
^10+ 





dx^{x)'^/3k{n)fk{x) 

k=l 
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With this provisions in mind, it is then only a question of sheer calculation (see 
Appendix C) to show that the pm^n(y,z) for 2 : > 0 as given in (16), (17), (18) 
and (19) correctly converge to the values (37), (38) and (39) for every possible 
ordering of n, m. For instance for n > m > 0, with x = aw and recalling also 
that 7rfc(0) = 6kfl (so that 7r„_m(0) = 0 because n > m), in the limit O’*" we 
immediately have 


Pm,n(|/,0^ 


(>0+ 


= '^m{y) / dx 
Jo 


and so on for the other two cases 


a" - - (a" - a"+i) 7rfc(0) 

n—m 

(a- - a-+i) 4,( 


k=0 


6{x) 




k=0 


= 0 


□ 


Proposition 5.4. The terms Q, A, B and C in the Proposition 5.3 can be expressed 
in terms of finite combinations of elementary functions: when z < 0 (namely aps > 
\t) we have for m > n >0 

m m , ^ \j—k / ‘\ ^ 

Qm,niy^^) = - ~j - 

k=n j=k ' 7 £—Q 

Tim-jiy - z)Tij+fiay) <h(j + 1; j + 7 + 1; ay) 

Here and in the following <F(q;; /3;x) are confluent hypergeometric functions. When 
instead z > 0 (namely aps < Xt) we have for every n,m >0 


Am,n{y 1 z) T^miy ^) ^ ^ fikijd) 

while for n > m > 0 it is 


k=0 


1 + TTkiaz) - ^ 'n-jiaz) 
j=o 


(26) 


n—m n+1 


Bm,n{y.z) = ^ TTfc (^) + 1) 


(27) 


A:=0 


1=0 


[az 


Yk\ 


<h (£, /c + £ + 1 , (1 — a)z) 


{k + i)\ 

and for m > n > 1 (for m > n = 0 we always have Cm,o{y) z) = b) it is 


„—{l—a){y—z) ^ 

Cm,n{y-, Z) = — ^ fifin) ^ 

(.=1 k=n j=0 


-i-ifk + i-j - 1 


7rj(ay)7rm+e-jia(y - z))^{k + i - j,m + i - j + l,a{y - z)) 


(28) 



N Cufaro Petroni and P Sabino: Correlated Poisson processes 


19 


Finally, since the parameters a, (3 of the $( 0 ;, x) involved in the previous equations 
are integer numbers with 0 < a < (3, our confluent hypergeometric functions are just 
finite combinations of powers and exponentials according to the following formulas 


$(0,/3,a;) = 1 


7r^-i(a^) 

TTg_i{x) 


Ida—7+1 ( 2 ^) 


f3 > a = 0 
fl > a > 1 


Proof: The detailed proof unfolds along a sequence of integrations based on known 
results and is here omitted for the sake of brevity (see Appendix D) □ 

We end this section with a short list of a few explicit examples of joint probabilities 
holding in the region afis > Xt: 

Po,o{s,t) = 

Pifi{s,t) = -[(1 - a)(l - e~^*) + a/as - Xt] 

ct 

p — flS 

Pipi-sfl) = -[At - (1 - a)(l - 

Qj 

P 2 ,o(s, t) = [ 2(1 - a)(l + ap,s){l - + (a/is - At)^ - 2(1 - a)Xt] 

P 2 U-s, t) = —— [(1 — a)(a — 4 — (1 — a)At — a/is)(l — 

’ 

+ At(a^ — 5a + 4 + ap,s — At)] 

g 

P 2 , 2 {s,t) = ^^[ 2(1 - a)(3 - a + (1 - a)At)(l - e“^‘) 

+ At(At-2(l -a)(3-a))] 


6 Cross-correlations and relative timing 

In this section we will briefly discuss the main correlation properties of the two 
processes. We hrst of all look at the point processes and we remark that, by recasting 
the self-decomposability equation (9) in the form 

Xk = ^Wk + Bfll)Zk Wfc = - W ~ (Biiia) 

X ja 

the point processes appear as 

n n 

k=f) k=0 

where now ~ ^i(A) and Wk ~ Si(h) play at once the role of the correlated 
renewals. It is interesting to point out then that, at variance with other models [6], 



N Cufaro Petroni and P Sabino: Correlated Poisson processes 


20 


- s=0 s=30 - - 3=70 s=110 - s=140 



Figure 6: Cross-correlation p{s,t) of the two Poisson processes M{s) and N{t) 
estimated by Monte Carlo simulations. Here a = 0.5, and A = p = 20, namely 
= 0.5 < 1 


we are no longer tied to take truly coincident shocks: we will show indeed that 
with non-zero probabilities the values of the paired, and correlated renewals X^, H4 
(waiting times) can be in an order whatsoever, and they would almost never coincide. 
As a consequence the propagation of the shocks from a process to the other will 
quite plausibly happen with delays whose random sizes (and directions) could also 
be modeled by suitably choosing our parameters a, A and p. And moreover the 
random times and S'„ will be correlated by the summing up of the renewals, but 
will never fall at the same instant. This relative timing apparently allows for an 
enhanced flexibility of the model in the practical applications because we no longer 
have to rely on common shocks, but rather on correlated and randomly delayed ones 


More precisely we can now single out two possible regimes for our processes: 
“^/a < 1 and “^/a > 1. It is then easy to see that for every fc = 1, 2,... 


Xk 


a/i 

T 


Wk + Bk{l)Zk>-^Wk>Wk 



1 


and hence we first of all have 


P{Xu>Wu] = l -f>l 


On the other hand for “^/a < 1 the probability P {Xk > H4} can still be explicitly 
calculated by taking into account the laws specified in the Section 4, and in this 
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case it is possible to show that 


P{Xk>Wk} 


(1 — a)fi 
A -|- (1 — o)yU 


a/i ^ 

T - 


1 


a value ranging from 0 to 1 according to the different possible choices of the param¬ 
eters a, A and /i 

As for the relative timings T„, Sm of the shocks along the point processes them¬ 
selves, an explicit calculation of P {Tn < Sm} is certainly possible, but its results 
would turn out to be rather cumbersome because it would involve two or three con¬ 
volutions of (positive and negative) Erlang laws with different parameters. We will 
then conhne ourselves here to produce just the cross-correlations between Tn, Sm- 
since it is easy to check that 


cov[Xk,W,] = ^6ke 

An 

it is also apparent that 

n m n m 

COV [Tn, Sm] = EE COV [Xk, We] = — EE 

k=l e=i ^ k=l i=l ^ 

and hence the cross-correlation coefficient of T„, Sm will simply be 

m 
m 

Finally, even closed formulas for the cross-correlation coefficient p{s,t) between 
the Poisson processes M{s) and N{t) are still derivable on the ground of our previous 
results about the joint distributions, but it would be too long to thoroughly elaborate 
them here. As an alternative we have chosen to show the plots of their estimates 
based on a sample of 10® Monte Carlo simulations of their trajectories as shown 
in the Figures 6 and 7. There the behavior is displayed of p{s,t) as a function 
of t for different, hxed values of s. More precisely, in the Figure 6 we have taken 
a = 0.5 and A = /x = 20 as the values for the relevant parameters of our coupled 
processes (then we have “^/a < 1), while in the Figure 7 the same parameters are 
a = 0.8, A = 20 and /i = 40 (and then “^/a > 1). It is apparent from these pictures 
that the behavior of p{s, t) is comparable to that of the self-correlation of a simple 
Poisson process, but for the fact that the cumulative effect of the correlate renewals 
produces a smoothing of the shape around the maximum values near t = s. At first 
sight this could look as a little difference, but in the domain, for instance, of the 
hnancial applications even small deviations among the models could produce huge 
differences in gains and losses 


n Am 


^ n.m. ^ 


a/V 

av^ 


m for n < 
for n > 
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- s=0 s=30 - - 3=70 s=110 - s=140 



Figure 7: Cross-correlation p(s,t) of the two Poisson processes M(s) and N(t) 
estimated by Monte Carlo simulations. Here a = 0.8, A = 20 and p = 40, namely 
= 1.6 > 1 

7 Conclusions and further inquiries 

It is apparent that, within the model discussed in the Section 2, from the self- 
decomposability alone we can only get pairs of rv’s X^Y with positive correlations 
0 < rxY < 1 steered by the value of a parameter a. It would be interesting however 
to widen the scope of our models in order to achieve also Poisson processes whose 
correlation can span over all its possible values (both positive and negative) by 
changing the value of some numerical parameter. In this respect it is important to 
remark - as pointed out in the Section 3.3 - that while two rv's X and Y which 
are, for instance, marginally exponentials can also be totally correlated {rxY = 1), 
they can not instead be totally anti-correlated {rxY = —1) because this would imply 
some linear dependence with a negative proportionality coefficient, and that would 
be at odds with the fact that both our rv^s take arbitrary large, but only positive 
values. Hence two exponential rv’s X and Y can always have a negative correlation, 
but only up to a minimal value which in any case must be larger than —1. We also 
showed in the Section 3.3 that this minimum is reached when between X and Y 
there is a peculiar kind of mutual functional, decreasing dependence, albeit clearly 
not a linear one. A model to produce pairs X,Y of rv’s which are marginally 
exponentials, and which - following the value of a numerical parameter a - show all 
the possible correlation values will be discussed in a subsequent paper 

Our results in any case show that the self-decomposability, joined with the sub- 
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ordination techniques, can be a promising tool to study dependency beyond the 
Gaussian-Itb world. We have shown indeed how to obtain dependent exponen¬ 
tial (gamma) rv^s that can be used to create and simulate dependent Poisson pro¬ 
cesses without resorting to dehnitely coincident jumps, but the path is now open 
to produce more general dependent gamma (Erlang at hrst) rv^s to simulate de¬ 
pendent variance gamma processes. A further extension could then be to study the 
self-decomposability of density functions that have a strictly proper rational char¬ 
acteristic function (Exponential Polynomial Trigonometric, EPT laws) in order to 
construct 2-dimensional correlated EPT rv’s (see for instance [16]). Finally it would 
be expedient to explore the Markov properties of the 2-component Poisson processes 
(M, {t),N{t)) with dependent marginals that we have introduced in this paper and 
the Master equations ruling them: this too will be the subject of future inquiries 
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A Notations 


All along this paper we will adopt the following notations: for a Poisson law ip(Q;) 
we will introduce the symbols 


^n(®) T 6 

n\ 


nn(a) = 'Y^'Kkia) 


a > 0 n = 0,1, 2 


k=n 


and for a binomial law f8(n, 1 — a) the notation 


f3k{n) = 


n = 0, 


k = 0 


0 < a < 1 


, n 


(-)an-fc(l_a)fc n = l,2,..., k = 0,..., 

It will be understood moreover that 

7rn(0+) = 5n,0 

We will also use bot the Heaviside function d, and the Heaviside symbol 0 


d{x) = 


1 a; > 0 

0 a; < 0 


0,= 


1 J >1 

0 j <0 


j = 0,±l,±2,... 


The pdf and chf of a standard Erlang law (£n(l) moreover will be denoted as 

(5(a;) 


fnix) = 




Tn{u) = 


1 

1 


(t^) 


n = 0 

n = 1,2,... 
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where it is understood for the Dirac delta d(x) that for every 6 > 0 


S(x) dx = 1 


Remark also that apparently 


lim / S(x) dx = 0 
2 ^ 0 + 


fn(x) = 7rn-i(x)d(x) n = 1, 2,... 

We will dually dedne for later convenience the functions 

n 

hn{x) = '^ldk{n)fk{x) n = 0,l,2,... 
k=0 

which are the pdf’s of the mixtures of Erlang laws €^(1) with binomial 23(n, 1 — a) 
weights for their indices k 


B A proof of Proposition 5.1 

To evaluate P {S'm < p, Sn < t} we drst remark that 

P {Sm < P: Sn < t} = P {M(p) > m, M(t) > u} 

and then that, being a Poisson process, M{t) is non-decreasing: as a consequence 

m <n and r < p M{t) < M{p) hence {M(r) > n} C {M{p) > m} 

n < m and p < r M{p) < M{t) hence {M{p) > m] C {M(r) > n] 

In the case m < n we then have for r < p 

P{M(p) > m, M(r) >n} = P {M(r) > n} 

while for p < t from the general properties of a Poisson process we get 

P {M(p) > m, M(r) > n} 

OO 

= ^ P {M(p) > m, M(r) > n \ M{p) = k} P {M(p) = k} 

k=m 

n 

= ^ P {M(r) > n I M(p) = k} P {M(p) = k} + P {M(p) > n} 

k=m 

n 

= ^ P {M(r - p)>n-k}P {M(p) = k} + P {M(p) > n} 

k=m 

In the same vein, when n < m we have for p <t 

P {M(p) > m, M{t) > n} = P {M(p) > n} 
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while for t < p we get 

P {M{p) > m, M{t) > n} 

OO 

= ^ P {M{p) > m, M{t) > n \ M{t) = k} P {M(r) = k} 

k=m 

m 

= ^ P {M(p) > m I M(r) = k} P {M(r) = k} + P {M(r) > m} 

k=n 

m 

= ^ P {M{p -T)>m-k}P {M(r) = k} + P {M(r) > m} 


k=n 


Remark that for m = n both the cases lead to the same result, namely 

pm(p)>n M(T)>n} = l whenr<p 

^ ^ ^ I -P {^(p) — when p < T 


that can also be conveniently summarized as 

P {M(p) > n, M(r) > n} = P {M(p At) >n} 

On the other hand for m <n we have 
for r < p P {M(r) > n} 

n—1 

for p < r P {M(p) >n} + ^ P {M(r — p) >n — k}P {M(p) = k} 

k=m 

that can also be put in the form 

n—1 

P {M(p A r) > n} + P(r - p) ^ P {M(r -p)>n-k}P {M(p) = P} 

k=m 

while for m > n it is 


m—1 


for T < p 

P {M(r) > m} + ^ P {M(p - t) > m - k} P {M(r) = k} 

for p < T 

k=Ti 

P {M(p) > m} 

namely 



m—1 

P {M{p At) > m} + -d(p — r) ^ P {M(p — r) > m — /c} P {M(r) = k} 

k=n 

In both cases the hrst terms can expressed as P {M(p At) > m V n}, and in this 
form they also coincide with the previous result for m = n. On the other hand the 
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extra term with the sum (which is absent for m = n) must be taken in consideration 
either when we have both m < n and p < r, or when it is m > n and r < p. All 
these provisions can then be comprehensively taken into account in the formula 

P{Sm < P, Sn< r} 

= P {M{p A r) > m V n} + [Qn-m'd{r - p) + Qm-n'd{p - r)] ■ 

(mVn) —1 

P W{\p - t\) > {my n) - k} P {M{p ^T) = k} 

k=mAn 

which hnally takes the form of Proposition 5.1 by using the notations adopted in 
the Appendix A for the Poisson distributions and the Heaviside symbols 


C Proof details for Proposition 5.3 

C.l The case afis > Xt, namely z < 0 

We begin with the case aps > \t, namely 2 : < 0, by recalling also that for the inte¬ 
gration variable in (13) it is 0 < tc < y. As a consequence, when the Proposition 5.1 
in used in (13), we always have 


0 < r 


y — w y — z 
-<- = s = p 

p p 


(29) 


On the other hand, since the conditions of the Lemma 5.2 are met, we can also 
restrict ourselves to evaluate Pm,n{s,t) for 0 < n < m 


C.1.1 m = n > 0 

In this case from (13), (29) and from the Proposition 5.1 we hnd 


y-w 


qn,n{y,z) = a dwhn{aw)P<M 

Jo [ \ P 

ry 

qn+i,n{y,z) = a dwhn{aw) 

Jo 


> n 


P<!M( -—- 1 >n + l 


P<^M 


ry 


qn,n+i{y,z) = a dwhn+i{aw)P <M 

Jo { \ P 

ry 


p 

w — z 

p 

y-w 


>l\PlM 


y-w 

p 


= n 


qn+i,n+i{y,z) = a dwhn+i{aw)P<M 

Jo I \ P 


y-w 


> n + 1 


> n + 1 
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and then from (12) 


VnAv^z) 


a / dw hn{aw)P < M 


w — z 




= b\P{M 


a / dw hn{aw) 7ro{w — z) TTniv — w) 


(") 



( 30 ) 
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C.1.2 m > n > 0 

In order to make use of Proposition 5.1 we remark that now 

m>n+l m+l>n m+l>n+l 

so that from (13) and (29) 


qm,n{y^^) = a / dwhn{aw) 

Jo 

m—1 ^ 

k=n 

ry 

qm+i,n{y,z) = a dwhn{aw 


P<M 


y-w 

y 


> m 


w — z 




PlM 


> m — k \ P < M 


y-w 


y-w 

y 


= k 




> m + 1 


w — z 




PlM 


m X 

+ Y.p\m 

k=n ^ 

ry 

qm,n+i{y,z) = a dwhn+i{aw) 

Jo 

m—1 ^ 

+ ©m-n-l P < M 

k=n-\-l ^ 

ry 

qm+i,n+i{y,z) = a dwhn+i{aw) 

Jo 

U —^11 ^ 


>m + l — k>P<M 


y-w 

fi 


y-w 


y 

w — z 


> m 


>m — k>P<M 


PlM 


y 

y-w 

y 


y-w 

y 


> m + 1 


k=n+l 


W — Z 


> m + 1 — k> P I M 


y-w 

y 
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and hence from (12) (in a slightly simplihed notation) 


Pm,n = a dw hr, 

Jo 


P<!M ( -—- I =m 

h 




h 


m—1 ^ 

^p{m 

k=n ^ 


h 


w — z \ , \ ^ \ ^ ^ ( V — w i 

= m — k> P < M i - = k 


h 


h 


—a / dw hn+i 

Jo 


P<!M ( -—- I =m 

h 


h 


-P<M{^]>l\p\M(y^^=m 


Om—n—1 


m—1 ^ 

E 


k=n+l 


h 


=m-k\p\M(y^]=k 


h 


h 


ry 

= a dw 

Jo 


h 


{hn - hn+i)P { M ( 1 = 0 i P <! M f I = m 


h 


/i 


+hnP < M ( -—^ ]=m-n\plM( -—- \ = n 


m—1 ^ 

+ {hn — hn+l)0m-n-l M 

k=n+l ^ 


/i 


=m-k\p\M{y^\=k 


h 


h 


ry 

= a dw 

Jo 


h 


hn+iP <! M ( -—^ ] = m-n\plM ( ^—- \ = n 


m ^ 

+ {hn — hn+l) M 

k=n ^ 


/i 


w — z \ — It’, 

= m — A;>P<M| - - = k 


h 


h 


so that we hnally have in full notation 


Pm,n(l/, z) = a dw< hn+l{aw)7rra-n{u! - z)7Tn{y “ w) 


(31) 


+ [K(aw) - /!„+i(0Ul)] E’^m-‘(® “ - «>) 


k=n 


Remark that this formula correctly encompasses also the case m = n > 0. As stated 
in the Proposition 5.3, we can also concisely write the overall result in the form 


( \ = I Q^^ri{y,z) m = n>0 

Pm,n[y, Z) < z) 171 > fl > 0 


( 32 ) 
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where for m > n > 0 and < 0 we define 


y _^ 

Qm,n{y^ z) = a dw hn{aw) 22 - z)v:k{y - w) 


k=n 


C.1.3 A normalization check 

We provide here a quick check of the correct normalization of the joint probabilities 
Pm,n calculated in the previous section for aps > \t, namely z < 0 < y. First remark 
that with X = aw 

py ^ 

Qmfi{y,z) = a dwho{aw)'22'^m-k{w - z)7rk{y-w) 

k=0 

/ ay 

dx5{x)'22'^rn-k (^- - Z^ TTfc (?/“-)= X] (-^)7rfc (|/) 

fc =0 ^ ^ 


k=0 




) ^ y^i-z)^-’^ _ 

k\(m — k)\ m\ 

k=0 ^ ’ k=0 




\m—k 




m\ 


■{y - z)"' = TTmiy - z) 


Then, since here only the Pm,n with m > n > 0 do not vanish, we have (by neglecting 
the arguments y,z) 

m—1 


^ ^ Pm,n ^ ^ ^ ^ Pm,n Qo,0 “1“ ^ ^ 


m=0 n=0 


Qm^m “1“ E (Q m,n Qm, 


n+l> 


oo / m 


m=l 
m—1 


n=0 


^,71+1 


Qofi + EE Qm,n E Q m,' 

m=l \n=0 n=0 

oo / m m 

Qo,o + EE Qm,n ^ ^ Qm,n j — Qo,o + E« 


m,0 


m=l \n=0 
oo oo 

~ ^ ^ Qm^O ~ ^ ^ ~ 1 

m=0 m=0 

which conhrms the normalization 


n=l 


m=l 


C.2 The case a/is < Xt, namely y > z > 0 


In this case, since 0 < z < y and 0 < w < y, when we use Proposition 5.1 in (13) 
instead of (29) we hnd only 


0 < r 


y-w 

/i 


0 < p = s = 


y- z 

y 


(33) 


SO that p and r can now happen to be in an order whatsoever. As a consequence 
Lemma 5.2 does not hold, and we must consider all the possible choices of m,n 
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C.2.1 n > m > 0 

In order to make use of Proposition 5.1 we remark that now 


m < n + 1 


m+l<n m+l<n+l 


so that from (13) and (33) 


9m,n = a dw hn{aw)P < M 


y — z y — w 
- - A - - 


+a / dw hn{aw) 


n—1 ^ 

Y,p{m 


> n — k \ P < M 


9m+i,n = a dw hn{aw)P < M 


z — w 


y-z y-w\ 
- A - > n 


y-z 


+ aQn-m-i / dwhn{aw) 


n—1 ^ 

—^_L1 ^ 


k=m-\-l 


> n — k> P < M 


y-z 


9m,n+i =a dwhn+iiaw)P < M 


y-z y — w 
- - A - - 


> n + 1 


+ a dwhn+i{aw) 


n . 

Y,p{m 


9^+1,n+1 = a dwhn+i{aw)P < M 


>n + l — k>P<M 


y-z y — w 
- - A - - 


y-z 


>n + l 


+ a dwhn+i{aw) 


n . 

E 


k=ra-\-l 


Z — W 


>n + l — k>P<M 


y-z 
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and hence from (12) 


Pm,n = a dw hn{aw) P < M 


> n — m> P < M 


-a / dw hn+i{aw)P < M 


z — w 


y-z 


> n + 1 — m\P{M 


= PlM 


y-z 


= m> a dw hn{aw) I 1 — P < M 


z — w 


y-z 


< n — m 


-hn+iiaw) 1 - P < M 


z — w 


< n — m 


= PlM 


y-z 


= m> a dw < hn{aw)P < M 


z — w 


+ [hn{aw) - hn+i{aw)] - P |m ^ 


z — w 


= n — m 


< n — m 


which, by plugging in the explicit Poisson probabilities, can be written as 


Pm,n{y^ - z) a dw [hn{aw) - hn+i{aw)' 


+T^ 7 n{y-z)a / dw\hn{aw)7rn-miz - w) 


n—m 

-[hn{aw) - hn+i{aw)] ^ nk{z - w)| 


— d^m,n{yjZ') 2lm,,n+l(|/) “1“ Bm,n{yjZ') 1 (l/j 

where, with the notations (17) and (18) adopted in the Proposition 5.3, we have 
dehned 

Am,n{y,z) = 7lm{y-z)a[ dw hn{aw) 


Bm,n{y,z) = Tlmiy - Z) a dw hn+l{aw) TTkiz - W) 
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C.2.2 n = m > 0 

In this case from (13), (33) and from the Proposition 5.1 we hnd 
gn,n = a j^dwhn{aw)P |m A 

Qn+in = a f dw hn{aw)P \ M (- — -A- —>n + ll 

Jo I V h h / J 

+a J^dwhn{aw)pi^M > i|p|m = ^} 

qn,n+i = dw hn+i{aw)P |m >n + l| 

+a j^dwK+i{aw)pi^M =n 

qn+i,n+i = dw hn+i{aw)P |m ^ ^ a ^ >n + l| 


and then from (12) 

Pmn = a- f dw hn{aw)P \ M (- — -A -—^ ) = n\ 

Jo I V h h / J 

-aj^dwhniaw) (^1-p|m = o|^ P =n 

-aj^dwhn+iiaw)(^l~ P =o|^p|m^^^^ =n 

= a J dw [hn{aw) — hn+i{aw)] P |m ^^ = n| 

-a l^dwhniaw)P S^M =o|p|m(^^^^^ = n| 

-a j^dwhn+i{aw)P =o|p|m(^^^^^ = n| 


which with the Poisson probabilities become 


Pn,niy^ = '^niy - z) tt dw [hn{aw) - h„+i(aw)] 

Jo 

+'^n{y — z) a [ dw hn+i{aw)no{z — w) 


+a / dw hn{aw)TTo{w — z)'Kn{y — w) 

J z 

-^n^nidJi z) ^n,n+l(l/) z) + z') p Cn^nidJi z') 
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where, as in the Proposition 5.3, along with (17) and (18) we also introduced (19) 

py ^ 

Cm,n{y,z) = a / dw hnjaw) TTm-k{w - z)'Kk{y - w) 


C.2.3 m > n > 0 

In this case it is 


m > n + 1 


m+l>n m+l>n+l 


so that from (13) and (33) 


(im,n ~ dw hn{aw)P <^M ^ - A- - j >m 

ry 

+ a dw hn{aw) 


m—1 ^ 

^p\m 




Qm+in = a dw hn{aw)P < M ( -— -A- —>m + l 

Jo I V h h / 


+ a dwhn{aw) 


m ^ 

^p|m 




I I ( y — ^ y — 

Qm n+1 =a dw hn+i{aw)P < M ( -A - - > m 

Jo [ \ P P J 

ry 

+ Qrn-n-ia / dw hn+i{aw) 


m—1 ^ 

E 

— ^_L1 ^ 


^]>m-k\p Mr-^]=k 


k=n+l ^ \ y V \ 

f {y — ^ y — iv\ 

qm+i,n+i = a dw hn+i{aw)P < M -A - - > m + 1 

Jo I V h h / 


+ a dw hn+i{aw) 


m ^ 

E 

— ^_L1 V 


k=n-\-l 


^j>m + l-A;|p{M^^j=A; 
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and hence from (12) 


/ 7 n / \ 7 / \T T-, 1 7, ^ I y ~ ^ y ~ \ 

Pm,n = a aw [hnyaw) — nn+i[aw)\ P <M \ -A- I = m 


+a / dw hn{aw) 


m—1 ^ 

T.p{m 

^k=n ^ 


fj. fi 


w — z \ 7 \ ^ \ ^ f y — ui 

= m — k>P<M 


fi 


= k 


fi 


fi 




fi 


—a / dwhn^i{ctw) 


m—1 ^ 

Qm-n-l 

k=n-\-l ^ 


fi 


=m-k\plM^y-^ 


fi 


= k 


-P{M{ -—-] >i\pIm( -—- I = m 


fi j ) V \ y 

On the other hand, since for every t > 0 

P {M{t) > 1} = 1 - P {M{t) = 0} 

we can also write 

f (y — ^ y — ^ 

Pm,n = a dw [hn{aw) — h„+i(atc)] P < M { - -A - - I = m 


-a / dw [hn{aw) — hn+i{aw)] P < M 


p p 

y-w' 


fy m ^ T- 

+a / dw hn{aw)'^^P < M 

k=n ^ 

—a / dw hn+i{aw) P < M 


p 


= m 


w — z\ — It’ 

= m — k \ P { M 


p 

f w — z 


= k 


p 


= m — k >P < M 


p 

y-w 

p 


= k 


= a J dw [hn{aw) — hn+i{aw)] P ^ 


y- z 

p 


= m 


+a I dwl hn+i{aw)P M (-) = m — nl-P<|M^^ 


p 


= n 


p 


+ [hn{aw) - hn+i{aw)] 

' w — z 

p yvi ^ 

k=n 


m ^ 


p 


= m — k \ P I M 


y-w 

p 
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which when the Poisson probabilities are introduced goes into 

Pm,n(y, z) = TTmiy “ z)a / dw [K{aw) - hn+i{aw)] 

Jo 

ry 

+a / dw hn+i{aw)'!:,n-n{^ - z)'!:n{y - w) 


(36) 


ry _^ 

+a / dw [hn{aw) - hn+i{aw)] ^ TTm-k{w - z)7rk{y - w) 


k=n 




namely the result of Proposition 5.3 once the dehnitions (17) and (19) are taken 
into account 


C.3 Continuity in 2 ; = 0 

We will now compare the values of Pm,n{y, z) separately listed in the Proposition 5.3 
for z < 0 and z > 0, and we will show that they connect with continuity in = 0, 
in the sense that for every y > 0 

Pm,n{y-: 0 ) Pm,n{,yt 0 ) 


For z < 0 (namely afis > Xt) the results are given in the Lemma 5.2 and in formu¬ 
las (30) and (31), so that for 2 ; t 0“ we simply have 


Pm,nijj 1 0 ) 

Pn,n(l/,0") 

Pm,n{y-i 0 ) 


a / dw hn{aw) TTq{w) Tlniy — w) 

Jo 

a j dw\^n+l{aw)T[rn-n{w)T[n{y - W) 

m 

+ [hn{aw) — hn+i{aw)~\ E '^m—k {w)'Kk{y 


k=n 


n > m > 0 
n = m > 0 

m > n > 0 



(37) 

(38) 

(39) 


On the other hand, when z > 0 (namely afis < Xt) the results are given in formu¬ 
las (34), (35) and (36), where appears also as an integration limit so that it is 
advisable to premise a few formal remarks. If the integrand is a regular function 
^{x) we of course have with y > z > 0 


lim 

2 : 4 . 0 + 


^ {x) dx 


(>0+ 


^{x) dx = 0 


lim 

24,0+, 


^{x) dx = / ^{x)dx 


When instead the integrand contains the distribution S{x) = (5(a; —O’*") concentrated 
on a; = O'*', as indeed happens in every first term of hn{x) which is /3o{n)S{x) = a^5{x) 
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(see also Appendix A), we have 

/ z ^0+ py 

^{x)5{x)dx = J ^{0)5{x) dx = ^{0) lirn J ^{x)5{x)dx = 0 

As a consequence in the limit J, O’*", by neglecting the vanishing terms, the following 
integrals on [0, z] only retain the contribution of the term /c = 0 in 


lim / dx^(x)hn(x) 
z^0+ Jq 


r>0+ 


r>0+ 


dx^{x)^l3k{n)fk{x) = / dx^{x)/3o{n)fo{x) 


k=0 


^0+ 


a"' / dx ^{0)S{x) = a"^^(0) 


while on the other hand for the integrals on [z, y] we have 


ry 

lim / dx^(x)hn(x) 
zi0+ 


ry 

lo+ 


/ y 

dx^{x) E /3fc('n)/fc(x) / dx ^(yX^ ^ ^ (3kirb) fk^x^ 


k=0 


k=l 


where now the term /c = 0 is apparently missing 

Take hrst Pm,n{y, A for 2 ; > 0 and n > m > 0 as given in (34): from the previous 
remarks, with x = aw and recalling also that 7rfc(0) = Skfi (so that 7r„_m(0) = 0 
because n > m), in the limit J, O’*" we immediately have 


(>0+ 


Pm,n(|/,0+) = Tlmiy) / dx 


- a”+' + a^TTn-miO) - (a" - 0 "+') 7rfc(0) 

n—m 


k=0 


6{x) 


qU _ ^n+1 _ / n _ ^n+1 


k=0 


= 0 


which coincides with Pm,niy, 0“) in (37). Then, always with 2 ; > 0, consider the case 
n = m > 0 given in (35): now (since 7ro(0) = 1) in the limit 2 : O’*" we have that the 

hrst integral exactly compensates the term k = 0 missing in the sum of the second 
integral, so that 


Pn,n{yi 0 ) 


f>0+ 


T^niy) J dx [a'^ — + a'^’''Vo(0)] (5(a;) 

py ^ 

+a / dw7ro(w)7r„(|/- w) ^/3fc(n)/fc(aw) 
-^0 k=l 

py ^ 

aJ'Tin{y) + a / dw7ro{w)7rn{y - w)'^/3k{n)fk{aw) 

k=i 

fy n ^ 

a / dw7ro{w)7rn{y - w)'^f3k{n)fk{aw) 


k=0 


py 

a / dwTTo{w)TTn{y — w)hn{aw) 

Jo 
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which again coincides with Pn,n(y, 0 ~) in (38). Finally consider the case z > 0 and 
m > n > 0 which is given in (36): here with the same remarks as before in the limit 
O’*" we again have Pm,n(y, 0“) of (39) 


(-0+ 


Pm,n(l/,0+) = TTmiy) I dx S{x) 


+a J dw ^hn+liaw)TTrn-niw)TTniy - W) 

m 

+ [hn{aw) - hn+i{aw)]'^7rm-k{w)7rk{y - w) 

k=n 


= a”(l - a)TTmiy) +a dw{ hn+liaw)TTm-niw)TTniy - w) 

J 0+ 


k=n 


+ [hn{aw) - hn+i{aw)]'Y^'x^_k{w)'Xk{y - w) 

k- 

= al dw{ hn+l{aw)7lm-n{w)'n-n{y - W) 

+ [hn{aw) - hn+i{aw)]'^7rm-k{u!)7rk{y - w) 


k=n 


because it is easy to see that the hrst term exactly compensates the missing term 
in the second integral: 


ray 


dx 


(0)7r„(|/) + (a” - ^ Trm-k{0)T^k{y) 


k=n 

m 


S{x) 


f'ay 

dx a”(l - a) ^ dm-k,onk{y) = a”(l - a)nm{y) 


k=n 


D A proof of Proposition 5.4 

To hnd an explicit elementary formula for Qm,niy, z) (here 2 : < 0, namely aps > \t) 
let us begin by supposing m > n > 1: in this case - within our usual notations - 
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from (15) we have 


3-(j/-2) fy 


Qm,n(y^^) = - /^o(n)Maw)+ y^/3e(n)f£(aw) 


m / \ 

^(”)(u.-2r-'=(!,-u.)' 


g-(2/-2) r<^y 

'ITT'- a™ Jo 




dx a^S{x) + /3e{n) fe{x) 


.y2r]{x-azr-yay-xY 


m—k / \k 


n ( \m—K ( \ 

a {—az) (ay) 


^ /3,(n) n 

h(^-^y-Jo ' 


dxe-^x^-\x - az)^-yay - x)' 


On the other hand (see [17] 3.383.1) with v = ay — x it is 


dxe-^x^-\x - az)^-yay - x^ = 6““^ / dve^{ay - n)^-^[a(|/ - z) - v 


m—k^k 


r<^y _ u\ 

= e-^y j dv e^{ay - vy-^’^ ^ j “ z)]\-v) 


i( ^,\m—k—i 


= e-^y 


m — k 


m — k 




(—1)"^ ^[ci{y — z)y / dve^{ay — vY 




= . ]{-ir-^-ya{y-z)nayY^^-^ 

i=0 ^ * T 

B(£, m — i + 1) <h(m — i + l;m — i + i + 1; ay) 

where $(«; /3]x) is a confluent hypergeometric function (see [17] 9.2) and B(q;,/ 9) is 
the beta function (see [17] 8.38) with 


B(«,/3) = 


r(a)r(/3) (a-l)!(/3-l)! 

r(Q; + /5) (q; + /3 — 1)! 
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As a consequence we have 


Qm,n{,y 1 


m\ a' 


k=n 

n m—k 

£=1 i =0 


{—az)"^ ^ {ay) 


m — k 
i 


l)^-^-^[a{y - z)]\ay) 




{m — i)\ 
{m — i + t)\ 


^{m — i + — i + £ + l; ay) 


3-(y-d 


m\ 


k=n 


m—k 




e=o 


i =0 


m — k 


[a{y - z)Y{ay) 




- * + 1; m - i + £ + 1; ai/) 


{m — i + £)l 


™^ / 1 \m—k—i 


EE 

k=n i =0 


ly 


m — I 
k 


AW 


£=0 


'Ki{y - z)'Km-i+i{ay)^{m - i + l]m - i + £ + 1] ay) 
and with the change of the summation index j = m — i 


Qm,n{yi^) ^ ^ ^ ^ 




AW 


a^ \k, 

k=n j=k ^ ^ ^=0 

TTm-j{y - z)TTj+e{ay)^{j + 1; j + £ + 1; ay) 

as stated in the Proposition 5.4. We have proved this result by supposing m > n > 
1, but it is possible to check now by direct calculation that also its extension to 
m > n > 0 gives the right results for Qm,o with m > 0, and for Qn,n with n > 1. 
We have in fact from the dehnition (15) that for m > n = 0 it is 


Qofi{y-,z) = e / dw6{aw) = e = Tio{y - z) 

Jo 

ry £m\ 

Qm,o{y,z) = a I dw ho{aw)^ ( - z)"^ ^{y - w)'' 

k=o k / 


m\ 


'0 


p-iy-^)(y — z)^ py 

— a dw 5{aw) = TTmiy — z) 


m\ 


while for m = n > 1 we have with v = a{y — w) 


Qn,n{yi z) 


o-{y-z) ry 


n\ 

e-(y-2) 


a / dw 


aP5{aw) + E /Jk{n)fk{aw) 


n\ 


e~°-y 

{ay)^ +-y 

^ ' a^ ^ 


k=l 


A(A 


{y - wY 


^ay 


^ A-i)Wo 


P’v^iay - v)^-^dv 
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and since (see [17] 3.383.1) 

e"v"'{ay — v)^~^dv = n + l)<l>(n + 1; n + /c + 1; ay) 

we finally get 



Qn,n{]J)Z) = j (fly) 


n\ 


0-o.y 


= e 


y)Yl 


k=0 


a"- {n + k)\ 


n\{ayY 
{n + k)\ 

{ay)^^{n + 1; n + /c + 1; ay) 


1 + ^ I3k{n) + 1; n + /c + 1; ay) 

(X 

^=1 

(3k{n) n\ 


It would be easy to see now that all these results can also be derived as particular 
cases from our explicit expression (15) that can hence be considered as completely 
general 

We proceed now to calculate for 2 : > 0 (namely for \t > ays) the explicit 
elementary form of the terms A,B and C defined in (17), (18) and (19). With 
X = aw, and by taking into account [17] 8.350.1 and 8.352.1, from (17) we first have 
for m > 0 and n > 1 


Am,n{y,z) = 7rm{y-z)a dw 

Jo 


f3o{n)fo{aw) +f3k{n) fk{aw) 


k=l 


= T^miy - Z) 

= T^miy - z) laJ" + fJkin) 


^jkj^ r k-i -x^ 


k=l 


k-l 

1 -^Tijiaz) 
1=0 


= T^miv - z)'^f3k{n) 

k=0 


k 

1 + nk{az) - E 7rj{az) 

j=o 


Since however this result can be extended also to n - 

/ az 

5{x) dx 


0 giving the correct result 
= T^m{y - Z) 


Jo 

it can be definitely used for every possible value of n, m > 0 as recalled in (26) of 
the Proposition 5.4 
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Then, for n > m > 0 and x = aw, by using [17] 3.383.1 we have from (18) 


= '^miy-z)a dw hn+l{aw) TTkjz - w) 

k=0 

n—m 

,{y dxhn+i{ 

u—r\ ^ 0 


— '^n 


— TTn 


= TTr 


k=0 
n—m n+1 


X TTk Z - 


(d,,{n + l) / dxfi{ 

I —n e—(\ "'O 


fc=0 £=0 
n—m 


X 

X]Tlk\Z - 

a 


i{y -^)Y1 


k=0 

e 


a”+Vfc (z) 


/3,(n + 1) 


E fJiya - 
(P- 


= TT, 


akk\Z^^ J, 

n—m 

Ay-z) '^Tik{z) 1 -”+^ 


dxx^ ^{az — xye' 


k=0 

n+1 


I3e{n + 1 ) 7 ^^^ 4> (£, /c + £ + 1 , (1 - a)z) 


£=1 

n—m 


{k + t)\ 


n+1 


= TT, 


i{y X] + 1) 


[az 


Yk\ 


k=0 


£=0 


ik + i)\ 


^ {i,k + i + 1, {1 — a)z) 


because $(0, k + l,x) = 1. This result also coincides with (27) in the Proposition 5.4 
Finally, for m > n > 0 and 0 < z < y, from (19) we hrst remark that Cm,o{y, z) = 
0 because ho{aw) = 6{aw) is now peaked outside the integration interval [z, y]. Then 
for m > n > 1 (the term of the sum with i = 0 again vanishes for the same reason 
as before) with v = a{y — w), and u = a{y — z) for short, from (19) we have 


Cm,n{y^^) = / dvhn{ay - v)Y'^rn-k 


k=n 


u — V \ fv 

- TTfc - 

a J \a 


u — V 


PU _ _ 

Y^ein) / dv f(,{ay-v)Y^rn-k 

£=1 k=n 

n pu ™ 

Y^An) / dvfi{ay-v)Y 

i=l k=n 


o-iy-z) 


V 

TTfc ( - 
a 

m—k„,k 


{u — v)^ 


a'"’ 

a^m\ 


{m — k)\k\ 




1=1 


k=n 



m 

k 


,-(y-z)p-ay 


dv 




{ay - V 


X-l -ay+v 


-(M-t;)™-V 


a™m! 


(£- 1 )! ^ \k 

e=l ^ ' k=n ^ 



{^-l)\ 

1 

dv e^ay - v)^-\u - 
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and since it is (see [17] 3.383.1) 


dv e^(ay - vy-\u - 
£-1 


j=0 

£-1 

j=0 

£-1 

j=0 


J 




m—k 


e-i-j 


iayy(-l) 


B(m — k + l,k + £ — ^^{k + i — j,m + i — j + l,u) 


j\{i - 1 - j)\{m + i - j)\ 


ly 


{ayyu^^^ ^^{k + i — j,m + i — j + l,u) 


we can write 


Cm^ni^y 1 '^) 


m £—1 

£=\ k=n j=0 


k P i — j — 1 
k 


j\{mP£- j)\ 

^-(l-a)(y-z) ^ ^ ^-1 

£=\ k=n j=0 


^{k + i — j,m + i — j + l,u) 
e-i-jfk + i-3 - V 


k 


7rj{ay)7rm+i-j{a{y - z))^{k + i - j,m + i - j + 1, a{y - z)) 


as stated in (28) of the Proposition 5.4 

As for the confluent hypergeometric functions appearing in all these results they 
in fact boil down to hnite combinations of elementary functions as stated in the hnal 
formulas of the Proposition 5.4. The case n = 0 immediately follows indeed from 
the dehnitions (see [17] 9.210.1). On the other hand, when 1 < a < /d, with the 
change of variables y = x — z and by taking sequentially into account [17] 3.383.1, 
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8.350.1 and 8.352.1 we have 


$(q;,/ 9 + l,a;) = 


B(/3 - q; + 1, a)x^ Jq 
I3\e^ I 


(13 - a)\{a - l)\xf^ Jo 

/j!e^ ^ 

{m — a)\{a — 1) \ x^ 7 


e^z'^-Hx- zf-'^dz 


e y{x-yY ^dy 


q; — 1 


7|'_1'|0-7-l / 7//3-7-lp- 


a;^(— 1 ) 


y^ ^ e '^dy 


= e 


= e 


(_ 1 )« 7 


Q:—7—1 




OL—l 

^ (/3-a)!(a-7-1)! 7r/3(a;) 


as stated in the proposition 
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